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Abstract 

Consider non-linear time-fractional stochastic heat type equations of 
the following type, 

dtUt[x) = + ll~^[Xcr{u) F {t,x)] 

in {d + 1) dimensions, where v > Q, j3 £ (0,1), a £ (0,2). The operator 
clf is the Caputo fractional derivative while —( —A)“/^ is the generator 
of an isotropic stable process and is the fractional integral operator. 

The forcing noise denoted by F {t, x) is a Gaussian noise. And the mul¬ 
tiplicative non-linearity cr : R —>■ R is assumed to be globally Lipschitz 
continuous. 

Under suitable conditions on the initial function, we study the asymp¬ 
totic behaviour of the solution with respect to time and the parameter A. 

In particular, our results are significant extensions of those in [10], [24], 

[25] and [11[. Along the way, we prove a number of interesting properties 
about the deterministic counterpart of the equation. 

Keywords: Space-time-fractional stochastic partial differential equations; 
fractional Duhamel’s principle; Caputo derivatives; noise excitability. 


1 Introduction and main results. 

1.1 Background material. 

Recently, there has been an increased interest in fractional calculus. This is 
because, time fractional operators are proving to be very useful for modelling 
purposes. For example, while the classical heat equation dtUt{x) = Aut{x), used 
for modelling heat diffusion in homogeneous media, the fractional heat equation 
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d^Ut{x) = Aut{x) are used to describe heat propagation in inhomogeneous me¬ 
dia. It is known that as opposed to the classical heat equation, this equation is 
known to exhibit sub diffusive behaviour and are related with anomalous diffu¬ 
sions or diffusions in non-homogeneous media, with random fractal structures; 
see, for instance, [22]. The main aim of this paper is study a class of stochas¬ 
tic fractional heat equations. In particular, it will become clear how this sub 
diffusive feature affects other properties of the solution. 

Stochastic partial differential equations (SPDE) have been studied in math¬ 
ematics, and in many disciplines that include statistical mechanics, theoretical 
physics, theoretical neuroscience, theory of complex chemical reactions, fluid dy¬ 
namics, hydrology, and mathematical finance; see, for example, Khoshnevisan 
[17] for an extensive list of references. The area of SPDEs is interesting to math¬ 
ematicians as it contains a lot of hard open problems. So far most of the work 
done on the stochastic heat equations have dealt with the usual time derivative, 
that is ,5 = 1. Its only recently that Mijena and Nane has introduced time 
fractional SPDEs in [24]. These types of time fractional stochastic heat type 
equations are attractive models that can be used to model phenomenon with 
random effects with thermal memory. In another paper [25] they have proved 
exponential growth of solutions of time fractional SPDEs-intermittency- under 
the assumption that the initial function is bounded from below. A related class 
of time-fractional SPDE was studied by Karczewska [15], Chen et al. [6], and 
Baeumer et al [1]. They have proved regularity of the solutions to the time- 
fractional parabolic type SPDEs using cylindrical Brownain motion in Banach 
spaces in the sense of [9]. For a comparison of the two approaches to SPDE’s 
see the paper by Dalang and Quer-Sardanyons [8]. 

A possible Physical explanation of time fractional SPDEs is given in [6]. 
The time-fractional SPDEs studied in this paper may arise naturally by consid¬ 
ering the heat equation in a material with thermal memory. 

Before we describe our equations with more care, we provide some heuristics. 
Consider the following fractional equation, 

d^ut{x) = -v{-A)°‘^^ut{x) 


with P S (0,1) and Sf is the Caputo fractional derivative which first appeared 
in [5] and is defined by 


d^ut{x) 


1 

r(i-/3) 



drUr{x) 


dr 

{t — r)P ' 


( 1 . 1 ) 


If uq(x) denotes the initial condition to the above equation, then the solution 
can be written as 


ut{x)= Gt{x - y)uo{y)dy. 

jRd 


Gt{x) is the time-fractional heat kernel, which we will analyse a bit more later. 
Let us now look at 


d^ut{x) = -i^{-A)°‘/^ut(x) + f{t, x), (1.2) 
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with the same initial condition uo{x) and f(t,x) is some nice function. We will 
make use of time fractional Duhamel’s principle [29, 31, 30] to get the 
correct version of (1.3). Using the fractional Duhamel principle, the solution to 
(1.2) is given by 

utix) = [ Gt{x - y)uoiy)dy + f f Gt-rix - y)d^~^f{r,y)dydr. 
jR-i Jo Jr<‘ 

We will remove the fractional derivative appearing in the second term of the 
above display. Define the fractional integral by 

It fit) ■= {t-T)~>-^f{T)dT. 

For every /3 £ (0,1), and g G L°°{R+) or g £ C'(R+) 

d?It ait) = ait). 

We consider the time fractional PDF with a force given by /(t, x) = ll~^g{t, x), 
then by the Duhamel’s principle, the mild solution to (1.2) will be given by 

Utix) = Gtix - y)uoiy)dy + / Gt-rix - y)gir,y)dydr. 

JR-i Jo JR‘^ 

The reader can consult [6] for more information. The first equation we will 
study in this paper is the following. 

d^ut{x) = -i/(-A)“/^ut(a;) + ll~t^[X(j[ut{x)) W it,x)], x £ R'^, (1-3) 


where the initial datum uq is a non-random measurable function. W it, x) is a 
space-time white noise with x £ R'^, ll~t^ and cr : R R is a globally Lipschitz 
function. A is a positive parameter called the “level of noise". We will make 
sense of the above equation using an idea Walsh [32]. In light of the above 
discussion, a solution ut to the above equation will in fact be a solution to the 
following integral equation. 

Utix) = iguo)tix) + X [ [ Gt-six-y)aiusiy))Widyds), (1.4) 

Jo Jr'^ 


where 


iguo)tix) 



y)uoiy)dy. 


We now fix the parameters a and /3. We will restrict /3 £ (0, 1). The 
dimension d is related with a and /3 via 


d < (2 A /3 ^)a. 


Note that when fi = 1, the equation reduces to the well known stochastic heat 
equation and the above restrict the problem to a one-dimensional one. This is 
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the so called curse of dimensionality explored in [12]. We will require the fol¬ 
lowing notion of "random-held" solution. We will need d < 2a while computing 
the norm of the heat kernel, while d < a is needed for an integrability 
condition needed for ensuring existence and uniqueness of the solution. 

Definition 1.1. A random held t > 0,x G is called a mild solution 

of (1.3) if 

1 . Ut{x) is jointly measurable in t > 0 and x G 

2. \/{t,x) G [0,oo) X J^^Gt-s{x - y)a{us{y))W{dyds) is well-dehned 

in by the Walsh-Dalang isometry this is the same as requiring 

sup supE|Mt(a;)p < oo! 

xGS.<‘ t>0 

3. The following holds in 

ut{x) = {guo)t{x) + X [ [ Gt-s{x - y)a{us{y))W{dyds). 

Jo Js.<‘ 

Next, we introduce the second class equation with space colored noise. 

dfUt{x) = -v{-A)°'^^ut{x) + ll~^[Xa{ut{x))F{t, x)], x G R"*. (1.5) 

The only difference with (1.3) is that the noise term is now colored in space. 
All the other conditions are the same. We now briehy describe the noise. 

F denotes the Gaussian colored noise satisfying the following property, 

E[F(t, x)F{s, y)] = (5o(t - s)f{x, y). 

This can be interpreted more formally as 


Gov 


(fidF, / tj^dF ) = 


dx / dy(j)s{x)ips{y)f{x -y), (1.6) 


where we use the notation f (/>dF to denote the wiener integral of (j) with respect 
to F, and the right-most integral converges absolutely. 

We will assume that the spatial correlation of the noise term is given by the 
following function for 7 < d, 


f{x,y) 


1 

\x-y\'^' 


Following Walsh [32], we define the mild solution of (1.5) as the predictable 
solution to the following integral equation 


ut{x) 


{guo)tix) -I- A 


/ /' 

jR'i Jo 


Gt-s{x 


y)a[us{:y))F{dsdy). 


(1.7) 
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As before, we will look at random field solution, which is defined by 1.7. We 
will also assume the following 

7 < a A d. 

That we should have 7 < d follows from an integrability condition about the 
correlation function. We need 7 < a which comes from an integrability condition 
needed for the existence and uniqueness of the solution. 

We now briefly give an outline of the paper. We adapt the methods of proofs 
of the results in [10] with many crucial nontrivial changes. We state main results 
in the next section. We give some preliminary results in section 3, we prove a 
number of interesting properties of the heat kernel of the time fractional heat 
type partial differential equations that are essential to the proof of our main 
results. The proofs of the results in the space-time white noise are given in 
Section 4. In Section 5, we prove the main results about the space colored noise 
equation, and the continuity of the solution to the time fractional SPDEs with 
space colored noise. Throughout the paper, we use the letter C or c with or 
without subscripts to denote a constant whose value is not important and may 
vary from places to places, li x € K'^, then |a;| will denote the euclidean norm 
of a; €, while when A C |A| will denote the Lebesgue measure of A. 

1.2 Statement of main results. 

Before stating our main results precisely, we describe some of the conditions we 
need. The first condition is required for the existence-uniquess result as well as 
the upper bound on the second moment of the solution. 

Assumption 1.2. • We assume that initial condition is a non-random bounded 

non-negative function uq : —>■ M. 

• We assume that cr : R —^ R is a globally Lipschitz function satisfying 
cr(x) < To-lx I with L„ being a positive number. 

The following condition is needed for lower bound on the second moment. 

Assumption 1.3. • We will assume that the initial function uq is non¬ 

negative on a set of positive measure. 

• The function a satisfies cr(x) > la-\x\ with l^- being a positive number. 

Mijena and Nane [24, Theorem 2] have essentially proved the next theorem. 

We give a new proof of this theorem in this paper. 

Theorem 1.4. Suppose that d < {2 A fi~^)a. Then under Assumption 1.2, 
there exists a unique random-field solution to (1.3) satisfying 

sup E|mj(x)|^ < forall t > 0. 

Here ci and C 2 are positive constants. 
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Remark 1.5. This theorem says that second moment grows at most expo¬ 
nentially. While this has been known [24], the novelty here is that we give a 
precise rate with respect to the parameter A. Theorem 1.4 implies that a ran¬ 
dom field solution exists when d < (2 A It follows from this theorem 

that TFSPDEs in the case of space-time white noise is that a random field so¬ 
lution exists in space dimension greater than 1 in some cases, in contrast to 
the parabolic stochastic heat type equations, the case /3 = 1. So in the case 
a = 2, /3 < 1/2, a random field solution exists when d = 1,2, 3. When /3 = 1 a 
random field solution exist only in spatial dimension d = 1. 

The next theorem shows that under some additional condition, the second 
moment will have exponential growth. This greatly extends results of [11], [13], 
[10] and [7]. 

Theorem 1.6. Suppose that the conditions of Theorem l.f are in force. Then 
under Assumption 1.3, there exists a T > 0, such that 

inf E|Mt(a;)P > forall t>T. 

Here C 3 and C 4 are positive constants. 

The lower bound in the previous theorem is completely new. Most of the 
results of these kinds have been derived from the renewal theoretic ideas devel¬ 
oped in [11] and [13]. The methods used in this article are completely different. 
In particular, we make use of a localisation argument together with heat kernel 
estimates for the time fractional diffusion equation. 

Remark 1.7. The two theorems above imply that, under some conditions, 
there exist some positive constants 05 and gq such that, 

a^) 2 o./{a- 0 d) < liminf ilogE|M 4 (x)p < limsup ^ logE[itt(a;)p < 

t^oo t t^ca t 


for any fixed x G 

The exponential growth of the second moment of the solution have been 
proved under the assumption that the initial function is bounded form below in 
[25]. This exponential growth property have been proved by [11] when /3 = 1 
and d = 1 when the initial function is also bounded from below. When /3 = 1, 
and the initial function satisfies the assumption 1.3, this was established by [10]. 
Chen [7] has established intermittency of the solution of (1.3) when d = 1 and 
f3 € (0,1) and ft G (1, 2) with measure-valued initial data. 

We will need the following definition which we borrow from [16]. Set 


St(X) := 


E|Mt(x )[2 dx. 
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and define the nonlinear excitation index by 


e{t) := 


lim 
A—>^00 


log loggt(A) 
log A 


The next theorem gives the rate of growth of the second moment with respect 
to the parameter A, which extends results in [10]. We note that for time t large 
enough, this follows from the theorem above. But for small t, we need to work 
a bit harder. 


Theorem 1.8. Fix t > 0 and x £ we then have 

loglogE|Mt( 3 :)p _ 2 a 
A^-oo log A a — dp 

Moreover, if the energy of the solution exists, then the excitation index, e(t) is 
also equal to 

Note that for the energy of the solution to exists, we need some assumption 
on the initial condition. One can always impose boundedness with compact 
support. 

The following theorem is essentially Theorem 2 in [24]. We only state it to 
compare the Holder exponent with the excitation index. This shows that the 
relationship mentioned in [10] holds for this equation as well: 77 < l/e{t). Hence 
showcasing a link between noise excitability and continuity of the solution. 

Theorem 1.9 ([24]). Let 77 < {a — Pd)/2a then for every x G {ut{x),t > 0}, 
the solution to (1.3) has Holder continuous trajectories with exponent 77 . 

All the above results were about the white noise driven equation. Our first 
result on space colored noise case reads as follows. 

Theorem 1.10. Under the Assumption 1.2, there exists a unique random field 
solution ut of (1.5) whose second moment satisfies 

sup E|ut(a:)p < C5 exp(c 6 A^“^^““'’'^^t) for all t > 0. 

a:eR‘^ 


Here the constants 05,05 are positive numbers. If we impose the further require¬ 
ment that Assumption 1.3 holds, then there exists a T > 0 such that 

inf E|Mt(a;)|2 > 07 exp(c8A2“/(“-T'^)t) for all t > T, 

xGB(0,t'3/“) 

where T and the constants 07,03 are positive numbers. 

Remark 1.11. Theorem 1.10 implies that there exist some positive constants 
og and oio such that 

cgA2«/0-/37) < hminf ilogE|Mt(a;)[2 < limsup - logE|Mt(x)[2 < 
t->oo t t^oo t 

for any fixed a: G 
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Theorem 1.12. Fix t > 0 and x G we then have 

loglogE|ut(a;)p 2 a 

iim - - --- = - 

A^oo log A a — 7/3 

Moreover, if the energy of the solution exists, then the excitation index, e(t) is 
also equal to ^ . 

We now give a relationship between the excitation index of (1.5) 

Theorem 1.13. Let rj < {a — fi'y)f2a then for every x G {ut{x),t > 0}, 
the solution to (1.5) has Holder continuous trajectories with exponent rj. 

While the general strategy of proof is the same as that used in [10], here we 
develop some new important tools. For example, we need analyse the heat kernel 
and prove some relevant estimates. In [10], this step was relatively straightfor¬ 
ward. But here the lack of semigroup property makes it that we need to work 
much harder. To address this, we heavily rely on subordination. This insight, 
absent in [7[ makes it that we are able to vastly generalise the results of that 
paper. Another key tool is showing that with time, {QuQ)t{x) decays at most 
like the inverse of a polynomial. This also requires techniques based on subor¬ 
dination. 


2 Preliminaries. 

As mentioned in the introduction, the behaviour heat kernel Gt{x) will play an 
important role. This section will mainly be devoted to estimates involving this 
quantity. We start by giving a stochastic representation of this kernel. Let Xt 
denote a symmetric a stable process with density function denoted by p{t, x). 
This is characterized through the Fourier transform which is given by 

p(L?) = e-‘'^lf“. (2.1) 

Let D = {Dr, ?■ > 0} denote a /3-stable subordinator and Et be its first 
passage time. It is known that the density of the time changed process XEt is 
given by the Gt{x). By conditioning, we have 

poo 

Gt{x)= p{s, x)fEt{s)ds, (2.2) 

^0 

where 

fstix) = (2.3) 

where g^(-) is the density function of Di. and is infinitely differentiable on the 
entire real line, with gp{u) = 0 for m < 0. Moreover, 

gp{u) ^ exp{—11 — /3|(it//3)^/^^“^^} as u —> 0-I-, (2.4) 





and 


5/3 (m) 


r(i-/3) 


,-/3-l 


as u —>■ oo. 


(2.5) 


While the above expressions will be very important, we will also need the 
Fourier transform of Gt{x). 


where 


and 


oo 

Ep{.x) = ^ ■ 


fe=0 


r(i + /3fc)’ 


( 2 . 6 ) 


i + r(Lwx ^ ^ i + rdVffl-'x ^ 


see, for example, [26, Theorem 4]. Even though, we will be mainly using the 
representation given by ( 2 . 2 ), we also have another explicit description of the 
heat kernel. 

Using the convention ~ to denote the Laplace transform and * the Fourier 
transform we get 


Gl{x) 


A/3 + 


( 2 . 8 ) 


Inverting the Laplace transform yields 


Gm = 


(2.9) 


where 


Ep{x) = ^ 


fc =0 


T{l + Pk)’ 


( 2 . 10 ) 


is the Mittag-Leffler function. In order to invert the Fourier transform, we will 
make use of the integral [14, eq. 12.9[ 


/*oo ^ 

/ cos{ks)E0^a{—as^)ds = 7^3 3 


(l.l),(a./9),(l,M/2)- 

/o ’ k ' 

Qj 



where TZ{a) > 0, /3 > 0, fc > 0, a > 0, is the H-function given in [19, 
Definition 1.9.1, p. 55[ and the formula 

1 /.oo ^ 

- = - / f{0cos{^x)d^. 

J — OO ^ J 0 

Then this gives the function as 


Gt{x) = 7 —j-ff 


2,1 

3.3 


>r 

(l,l).(l,/3).(l.a/2)- 

_ l^t/3 

(l,a).(l,l),(l,a/2)- 


( 2 . 11 ) 
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Note that for a = 2 using reduction formula for the H-function we have 


1 


Gtix) = 


xp 

(b/3)-| 

lytl^ 

(1.2)- 


Note that for /3 = 1 it reduces to the Gaussian density 

1 / UI2 


( 2 . 12 ) 


Gtix) = 


■ exp 


F, 

4z/t 


{Ai/TTtyG 

We will need following properties of the heat kernel of stable process. 


(2.13) 


pit, x) = t-'^/“p(l, t-i/“x). 

Pist, x) = s-'^/°‘pit, s-i/“x). 

• pit, x) >pit, ?/)whenever |x| < |?/|. 

• For t large enough so that pit, 0) < 1 and r > 2, we have 

pit, -(x - y)) > pit, x)pit, y). 

T 

All these properties, except the last one, are straightforward. They follow from 
scaling. We therefore provide a quick proof of the last inequality. Suppose that 
t is large enough so that pit, 0) < 1. Now, we have that < EM v EM < 
|x| V |?/|. Therefore by the monotonicity property of the heat kernel and the fact 
that time is large enough, we have 

pit, -ix-y)) >pit, |x| V |y|) 

T 

>pit, |x|) Apit, |?/|) 

>pit, \x\)pit, |?/|). 


We will need the lower bound described in the following lemma. The upper 
bound is given for the sake of completeness and is true under the additional 
assumption that a > d, a, condition which we will not need in this paper. 

Lemma 2.1. (a) There exists a positive eonstant ci such that for all x € 

(b) If we further suppose that a > d, then there exists a positive constant C 2 
such that for all x G 
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Proof. It is well known that the transition density p(t, x) of any strictly stable 
process is given by 


Cl 


\d+ct 


(2.14) 


<P((, A 

where ci and C 2 are positive constants. We have 

poo 

Gt{x)= / p(s, x)/i 5 ,(s)ds, 

which after using (2.3) and an appropriate substitution gives the following 

poo 

Gt{x) = / p{{t/u)^, x)g^{u)du. 

Jo 

Suppose that |x| < then t/\x\°‘^^ > 1. When we have u < we 

can write 


p{{t/u)^, x)gi 3 {u)du >C 5 {t/u)~^‘^^°‘gp{u)du 

Jo 

> ce f {t/u)~^'^^°‘gp{u)du (2-15) 

Jo 




u^‘^/‘^gp{u)du. 


Since the integral appearing in the right hand side of the above display is finite, 
we have Gtix) > whenever \x\ < We now look at the case 

|a;| > 



p{{tfu)^, x)gi 3 {u)du > 


Co-. 


{t/uY 


> 


> 




d-\-OL 


gp{u)du 


G f 


Cll 




(2.16) 


where we have used the fact that fY(u) ^gp{u)du is a positive finite constant 
to come up with the last line. 

We now use the fact that p{{t/u)^, x) < ci : we have 


Gt{x) < Cl 



y^Pd/a 

-[m^9p{u)du 


Cl 

iPdja 



U^'^l'^gp{u)du. 
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The inequality on the right hand side is bounded only \ia > d. This follows from 
the fact that for large u, gpiu) behaves like So we have Gt{x) < 

Similarly, we can use p{{t/u)^, x) < to write 


/■oo 


\x 


d-\-a. 


0 


u ^gp{u)Au. 


Since the integral appearing in the above display is finite, we have Gt{x) < 
°;b*, „ . We therefore have 




Gt{x)<C5(^t-^^/^ A 


tP 


Id+a 


□ 


Remark 2.2. When a < d, then the function Gt{x) is not well defined ev¬ 
erywhere. But its representation in terms of H functions, one can show that 
a: = 0 is the only point where it is undefined. We won’t use the pointwise upper 
bound. The lower bound is trivially true when a; = 0. 

The T^-norm of the heat kernel can be calculated as follows. 

Lemma 2.3. Suppose that d < 2q:, then 


[ GUx)dx = C*t-^‘^/‘^, 

aR<i 


(2.17) 


where the constant C* is given hy 

(z/)-‘'/“27r‘'/2 


G* = 


:r(|) (27r)« 


Q 


z‘^^‘^-\Ep{-z)fdz. 


Proof. Using Plancherel theorem and (2.9), we have 


\Gt{x)\'^dx = 


1 


\Gt{OM = 


1 


{27tY 

2Tr‘^G 1 

W Jo 

(z/t^)-‘^/“27r"*/2 1 

(27r)'' 


(27r)" 


/R'^ 


\Ep{-n\frt^)\^d^ 


«r(f) 


pOO 

/ r‘^-\E0{-iyr°‘tf^)fdr. (2.18) 

Jo 

o 

zP°^-\Ep{-z)fdz. (2.19) 


To finish the proof, we need to show that the integral on the right hand side of 
the above display is bounded. We use equation (2.7) to get 


ydj 0L—\ 


/o (i + r(i-/3)^)2 


dr < 


^d/a-l(^^(_^)) 2 d^ 
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■oo 


< 


ydj Ol — \ 


/o (i + r(i + /3)-i2)2 
Hence < oo if and only if d < 2a. 

Recall the Fourier transform of the heat kernel 

We will use this to prove the following. 

Lemma 2.4. For 7 < 2a, 


dz. 


( 2 . 20 ) 

□ 

( 2 . 21 ) 




( 2 . 22 ) 


where Cl = z'^/°‘ ^{Ep{-z))‘^dz. 

Proof. We have 




dC = 


2'k‘^I'^ 


de 


ill poo -I 

^ J r^-\Epi-i.r‘^t^)r^dr. (2.23) 


r(i) ./o 

(l/t^)-T'/“27r'^/2 


ar(f) (27r) ^0 


nOO 

, / z'^/^-\Ep{-z)rdz. 

do 


We used the integration in polar coordinates for radially symmetric function in 
the last equation above. Now using equation (2.7) we get 


7/a-i 


pOO 

L 


ylloL-l 


< 


/o (i + r(i + /3)-iz)' 


dz. 


(2.24) 


Hence z'>'/“ ^(i5^(—z))^dz < 00 if and only if 7 < 2a. In this case the upper 
bound in equation (2.24) is 


rl/a-l 


rdz = 


B( 7 /a, 2 - 7 /a) 


/o (l + r(l+/3)-iz)2““ r(l + /3)-T'/« ’ 
where B( 7 /a, 2 — 7 /a) is a Beta function. 

Remark 2.5. For 7 < 2a, 

B( 7 /a, 2 - 7 /a) ^ ^ ^ B{-f/a,2--f/a) 

rd-d)’/. -/ " (BM-0)d^< rd + fl-v= ■ 


□ 
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We have the following estimate which will be useful for establishing temporal 
continuity property of the solution of (1.5). 

Proposition 2.6. Let 7 < min{2,/3“^}a and h G (0, 1), we then have 

f dCds < 

Jo Jr-^ Kr ^ 

Proof. The computation in Lemma 2.4 we have 

JR'^ Is I 

- 2 [ Gt-s+hiOGt-AOj^d^ 

JR'* Is I 

jR-i l4r ^ 

Using integration in polar coordinates in R'^, and the fact that h{z) = 
Ep{—z) is decreasing (since it is completely monotonic, i.e. (—l)”h*^”^( 2 ;) > 0 
for all 2 > 0, n = 0,1, 2, 3, • • •), we get 

2 [ Gt-s+hiOGt-siOj^df 
jR<i Is I ' 

= 2 / Ep{-,.\^r{t-s + h)P)Ep{-,z\^r{t-sf)-^df 

jR<i Is I ' 

> 2 / Ef,i-u\^rit-s+hf)Ef,i-u\^rit-s+hf)j^d^ 

Jr-^ Is I 

= 2 G*(t-s + /i)-^^/“. 

Now integrating both sides wrt to s from 0 to t we get 

f [ |G,_,+;,(e)-G,_,(^)|2^_^d^dr (2.25) 

Jo Jr-^ Is I 

-GC(h)i-^T'/“ GJ‘l(t + h)i-^T'/“ Gi*U-^'>'/“ 

~ \ — f3"f/a 1—P^/a X—fJ^joL 

2 Gi*(t+ /i)i-^'^/“ 

I — Pjla l — (3d/a 

G*(/i)l-/37/a 

l — P^fa I — Pjla X—fJ^ja 

G*(/j)l-/37/a 

“ X — a ’ 

the last inequality follows since t <t!. 

□ 
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Lemma 2.7. Suppose that 7 < a, then there exists a constant ci such that for 
all X, y £ W^, we have 

( f Gt{x-w)Gt{y-z)f{z,w)dwdz<-^. 

Proof. We start by writing 

/ / p{t, X — w)p{t', y — z)f{z, w)dwd 2 : 

jRd jR<i 

= / p{t + t', X — y + w)\w\~'^ dw 

Jr'^ 


< 


C2 


(t + ' 

We use subordination again to write 


[ [ Gtix-w)Gtiy - z)f{z, w)dwdz 

Jr^ jRd 

p PCX) nOO 

= / / Pis,x-w)p{s',y-z)fxt{s)fxtis')dsds'fiz,w)dwdz 

JR'^xR'^io Jo 

POO PCiO p 

= / / pis, X-w)p{s', y - z)f{z, w)dwdzfxtis)fxtis')dsds' 

Jo Jo Jr’^xR'^ 

p(00 pCiO p 

= / / pis, X-w)pis', y - z)fiz, w)dwdzfxtis)fxtis')dsds' 

Jo Jo JR'ixR'^ 


^00 pOO 


< 


< 


C2 


:fxtis)fxtis')dsds' 


is + s')^/^ 

pCXD pCO 

/ / -^/^t(s)/^t(s')dsds'. 

Jo Jo ^ ' 


/o Jo 

^00 /•OO 


Recalling that fxtis') is a probability density, we can use a change of variable 
to see that the right hand side of the above display is bounded by 




Since the above integral is finite, the result is proved. 


□ 


The next result gives the behaviour of non-random term for the mild formu¬ 
lation for the solution. For notational convenience, we set 


{gu)tix):= Gtix - y)uoiy)dy. 
Jr‘^ 


The proof will strongly rely on the representation given by (2.2) and we will 
also need 

iGu)tix):= pit, x-y)uoiy)dy, 
jR'i 
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where p{t, x) is the heat kernel of the stable process. We will need the fact 
that for t large enough, we have {Qu)t{x) > for x € B(0, We will 

prove this fact and a bit more in the following. The proof heavily relies on the 
properties of pit, x) which we stated earlier in this section. 

Lemma 2.8. Then there exists a to > 0 large enough sueh that for all t > 0 

{Qu)t+to{x) > cit~‘^^°‘, whenever x € 5(0, 

where ci is a positive constant. More generally, there exists a positive constant 
K > 0 such that for s <t and t > tg, we have 

{Qu)s+to{x) > C 2 t~^, whenever x G 5(0, 


C 2 is some positive constant. 

Proof. We begin with the following observation about the heat kernel. Choose 
to large enough so that p{to, 0) < 1. We therefore have 

p{to, x-y) = p{to, 2{x - y)/2) 

> p{to, 2x)p{to, 2y) 

= ^P(^o/2“, x)p{to, 2y). 

This immediately gives 

{Qu)to {x) = / p{t, X - y)uoiy) dy 

> Cip{tol2°‘, x) [ p{to, 2y)uo{y)dy. 


We now use the semigroup property to obtain 

igu)t+to{x) = / p{t + to, X-y)uo{y)dy 

= / p(,t. x-y){gu)to{y)dy 

Jr'^ 

> C 2 p{t + to/2, x), (2.26) 

This inequality shows that for any fixed x, {gu)t+to{x) decays as t goes to 
infinity. It also shows that 

(i7u)t+tQ(a;) > cot~‘^^°‘, whenever |x| < t^^°‘. 

This follows from the fact that p{t + to/2, x) > if |x| < The more 

general statement of the lemma needs a bit more work. 

{gu)s+to{x) > C 2 p(s + to/2, x) 

to 


> C3 

> C3 


2s + to 
to 

2s + to 


d/c 


p{to, x) 


djc 


p{to, t 


l/a\ 
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Since we are interested in the case when s < t and t > to, the right hand side 
can be bounded as follows 

The second inequality in the statement of the lemma follows from the above. □ 


Lemma 2.9. There exists a to > 0 and a constant Ci such that for all t > to 
and all x G B(0, we have 

{gu)s+t{x) > ^ for all s <t. 

Proof. We start off by writing 


{Gu)t{x)= Gt{x - y)uoiy) dy 
JR'^ 

= / pis, x-y)fxtis)dsuoiy)dy 

jRd Jq 

nOO 

= / (Gu) six) fxAs)ds. 

Jo 

After the usual change of variable, we have 

poo 

iGu)tix)= igu)^t/u)i^ix)gi3iu)du, 

Jo 

which immediately gives 

igu)tix)> / (^u)(t/„)p(x)g;3(M)dM. 

Jo 

The above holds for any time t. In particular, we have 

igu)t+six)> / (^M)((t+s)/'u)'3(a;)5/3(M)du. 

Jo 

We now that note that x € B(0, so we have x € B(0, t^^°‘fu) and hence 

for t large enough and s < t,we have iGu)(^(^s+to)/u)f> ix) > (f)^^ by the previous 
lemma. Combining the above estimates, we have the result. □ 


Remark 2.10. The above is enough for the lower bound given in Theorem 1.6 
and the lower bound described in Theorem 1.10. But we need an analogous 
result for the the noise excitability result which hold for all t > 0. Fix t > 0 
such that pit, 0) < 1 whenever t > t. For any fixed t > 0, we choose k large 
enough so that 2^t > i. Set t* := 2^t and s = 2~^. 

pit, x-y) =pist*, x-y) 

= s-^/°^{x-y)) 

= s~‘^/°‘pit*, - 2y)). 
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For any fixed f > 0, we choose k large enough so that 2^t > i. 

p{t, x-y)> 2 s-i/“ 2 /) 

= 2‘^'=/X2'=t, 2i+'=/“x)p(2'=t, 2i+'=/“2/). 

Note that the above holds for any time t. We therefore have 

{Guo)ta+s{x) = / p{tQ + s,x-y)uo{y)dy 

> + s), 2i+'=/“x) [ p{2>^ito + s), 2^+’^/^y)uo{y) dy. 

Jr‘‘ 

We have that to + s > ^o- Therefore, 

p{2\to + s),2^+'^/^x)> ' P(2%, 2i+'=/“x) 

V s + ^0 y 


We thus have 

(euo)*„+.(x) > 2'^'=/^ 


^0 \ 


2dja 


■J 


+ to ^ 

So now since |x| < we have 


p{2Ho,2^+'^/^x) / p{2%,2^+’^/^y)uo{y)dy. 
Jr<‘ 


{Quo)to+s{x) > Cl 


1 


to + s 


2dja 


where the constant ci is dependent on to. We can now use similar ideas as in 
the proof of the previous result to conclude that if a; € B{0, we have 


iGuo)to+s{x) > C 2 


to + s 


2^dl a. 


Since we have s < t, we have essentially found a lower bound for (C/Mo)to+s(a;); 
a bound which depends only on t. This holds for any to > 0 and any t > 0. 

We end this section with a few results from [10]. These will be useful for the 
proofs of our main results. 

Lemma 2.11 (Lemma 2.3 in [10]). Let 0 < p < 1, then there exists a positive 
constant ci such that or all b > {ejp')^, 


E 

1=0 


5 \ J 

Tr 


> exp . 


Proposition 2.12. [Proposition 2.5 in [10]] Let p > 0 and suppose f{t) is a 
locally integrable function satisfying 


f{t) < Cl + K f [t — sY ^/(s)ds for all t > 

Jo 


0 , 
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where ci is some positive number. Then, we have 

f{t) < C 2 exp(c 3 (r(p))^/^K^/^t) for all t > 0, 
for some positive constants C 2 and C 3 . 

Also we give the following converse. 

Proposition 2.13 (Proposition 2.6 in [10]). Let p > 0 and suppose f{t) is 
nonnegative, locally integrable function satisfying 

f{t)>ci+K f {t — s)^~^ f{s)ds for all t > 0, 

Jo 

where ci is some positive number. Then, we have 

f{t) > C 2 exp(c 3 (r(p))^/^K^/^t) for all t > 0, 
for some positive constants C 2 and C 3 . 


3 Proofs for the white noise case. 


3.1 Proofs of Theorem 1.4. 

Proof. We first show the existence of a unique solution. This follows from a 
standard Picard iteration; see [32], so we just briefly spell out the main ideas. 
For more information, see [24[. Set 

:= iGuQ)t{x) 


and 

:= {guo)t{x) + X [ [ Gt-six - {y))W{dy ds) for n > 0. 

Jo JRd 

Define Dn{t ,x) := (x) — u[^\x)\'^ and Hn{t) := sup 3 ,gRd ,x). We 

will prove the result for t € [0, T], where T is some fixed number. We now use 
this notation together with Walsh’s isometry and the assumption on a to write 


D^{t,x) = X^[ [ Gl_,{x-y)¥.\<7{u^f^\y))^\y))\'^dyds 
Jo Jr’‘ 

<X^Llf Hn-i{s) [ G^t_,.{x - y) dy ds 
Jo Jr’‘ 


< X^L 


2 t2 


{t - syt^/^ 


ds 


We therefore have 


Hn{t) < X^Ll 



Hn-l{s) 

{t - syp/^^ 


ds. 
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We now note that the integral appearing on the right hand side of the above 

display is finite when d < a/fi. Hence, by Lemma 3.3 in Walsh [32], the series 
1 . 

Hn (t) converges uniformly on [0, Tj. Therefore, the sequence {u„} con¬ 
verges in and uniformly on [0, T] x and the limit satisfies (1.4). We can 
prove uniqueness in a similar way. We now turn to the proof of the exponential 
bound. From Walsh’s isometry, we have 

E\ut{x)f = \{guo)t{x)\'^ + [ [ Gt_^{x-y)E\a{usiy))\‘^dyds. 

Jo jR<i 

Since we are assuming that the initial condition is bounded, we have that 
< Cl and the second term is bounded by 

X^Ll j f G'^_,{x-y)E\us{y)\^dyds 
Jo Jr"^ 

r* \ 

< ciX^Ll / - -sup E|us( 2 /)pd?/ds. 

Jo H ~ 


We therefore have 

/■* 1 

sup E|Ms(a;)P < ci+C 2 X^LI / - -sup E|u^(i/)pds. 

a:eR‘i Jo U ~ ygR<i 

The renewal inequality in Proposition 2.12 with p = {a — df3)/a proves the 
result. □ 

3.2 Proof of Theorem 1.6. 

The proof of Theorem 1.6 will rely on the following observation. From Walsh 
isometry, we have 

E\ut{x)f = \{guo)t{x)\'^ + X"^ [ [ G‘f_^{x - y)E\a{us{y))\'^dyds. 

Jo Jr'^ 

For any fixed fg > 0, we use a change of variable and the fact that all the terms 
are non-negative to obtain 

E\ut+toix)\^>\{guo)t+toix)\‘^+>'^ll [ [ G^_,{x-y)E\us+toiy)\‘^dyds. 

Jo Jr"^ 

Using the above relation again, we obtain 

E|ut+t„(a;)p > \{guo)t+to{x)\^ 

+ X^ll [ [ G^_,{x-y)\{guo)s+to(.x)\^dyds 
Jo Jr-^ 

+ / / / / - y)Gl_^^{y - z)E\us^+toiz)\^dzdsidyds. 

Jo Jr'^ Jo Jr'^ 
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Using the same procedure recursively, we obtain 


E\ut+to{x)\^ 

> |(euo)t+io(a;)|^ 


2k ]2k 


k^l 




I 

Js.“ 


|(0'Uo)io+sfc(’2fc)| 


k 

— ^i) ^^k+l—i dSfc+1 —z, 

i=l 


(3.1) 


where we have used the convention that sq := t and zq := x. Let x € i3(0, 
and 0 < s < t and set 

(l?Mo)to+s(a;) > gt- (3.2) 

The existence of such a function gt is guaranteed by Lemma 2.9 and Remark 
2.10. We can now use the above representation to prove the following result. 


Proposition 3.1. Fix tg > 0 such that for t >0, 


E\ut+to{x)\^ > gt (A^Z^ci)' 




^ \ k{ty.—/3d) j a 


for xGB{0,t^/°), 


where ci is a positive constant. 

Proof. Our starting point is (3.1). Recall the notation introduced above, 


{guo)ta+s„{zk) > gt, 

whenever Zk € R(0,t^/“) and 0 < Sk < t. The infinite sum of the right of (3.1) 
is thus bounded below by 




fe=l 


0 Jr<‘ Jo 


/R<i 


/ TT -Si (zi-i,Zi) dzk+i-i dsfc+i_ 


We now make a reduce the temporal domain of integration and make an appro¬ 
priate change of variable to find a lower bound of the above display 




ntjk 


ntik 


nt/k 


k^l 


TT Gl. {zi-i,Zi) dzk+i-i dsfc+i_z. 

/o JM.'i Jo JB{0,tf>/°‘) 


We will reduce the domain of the function 

k 

Y[Gl.{zi-i,Zi), 


i=l 


by choosing the points Zi appropriately so that they are "not too far way". We 
choose zi G R(0, such that \zi — a;o| < In general, for i = 1, • • • , fc, 
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we choose Zi S fl i?(0, An immediate consequence of this 

restriction is that 




Cl 




2d^ /a ' 


Since the area of the set B{zi-i, sf^°‘) fl 5(0, t^/“) is C 2 sf^^“, we have 


d0fo 


nt / k p nt / k n nt / k n ^ 

/ / / / •■• / / ]^Gsi(^i-i)^i)dzfc+i-*ds 

Jo JR'i Jo JR<i JO JB(0, t/3/“) . , 


fc+1 —2 


> 


^t/fc ^t/fc fc 

I 


2=1 

,fc 

3 


'0 S, 


d/3/a 


dsi•••dsfe 


- c'= 

— C 3 


^ {a-dp)k/a 

k. 


Putting all the estimate together we have 


E|wt+t„(x)|2 > gl Y, ( - 


fc=0 


{a — d^)k/oi 


□ 

Proof of Theorem 1.6. We make the important observation that gt decays 
no faster than polynomial. After a simple substitution and the use of Lemma 
2.11, the theorem is proved. □ 

Remark 3.2. It should be noted that we do not need the full statement of 
Proposition 3.1. All that we need is the statement when time is large. 

3.3 Proof of Theorem 1.8. 

Proof. From the upper bound in Theorem 1.4, we have that for any x € 
E|Mt(x)p < forall t > 0, 


from which we have 


loglogE ut(a;) p 2a 

limsup-^-r- < - —. 

A^oo log A a- dp 


Next, we will establish a lower bound. Fix x S for any t > 0, we can always 
find a time Iq such that t = f — to + to and t — to > 0. If t is already large enough 
so that X G B{D,t^/°‘) then by Proposition 3.1 and Lemma 2.11 we get 


lim inf 

\—¥oo 


loglogE|Mt(a;)p 
log A 


> 


2a 

a — dp 
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Now \i X ^ we can choose a /t > 0 so that x € i?(0, 

Then we can use the ideas in Proposition 3.1 to end up with 

OO / •/- \ k{<^~^d)/OL 

E|ut+to(a:)|2 |'_j 

and the result follows from this using Lemma 2.11. □ 


4 Proofs for the colored noise case. 


4.1 Proof of upper bound in Theorem 1.10. 

Proof. The proof of existence and uniqueness is standard. For more information, 
see[32]. We set 

x) := iGuo)tix), 

and 

x) := {guo)t{x) + X f f Gt-s{x-y)a{u^’^\s, y))F(dyds), n>0. 
Jo Js.'i 

Define Dn{t,x) := x) — x)p, Hn{t) := sup,j.g][jd Dn{t,x) and 

T,(t,y,n) = y)) — y))|. We will prove the result for t € 

[0, T] where T is some fixed number. We now use this notation together with 
the covariance formula (1.6) and the assumption on a to write 


Drl(J'^ x) 

= /" / [ Gt-s{x-y)Gt-six - z)E[Y;{s,y,n)T,{s,z,n)]f{y,z)dydzds. 

Jo Js.'^ Jr’^ 

Now we estimate the expectation on the right hand side using Cauchy-Schwartz 
inequality. 

E[E(s,2/,n)S(s,z,n)] < T^E|u("')(s, y) 2 /)||m("')(s, z) z)\ 

/ \ 1/2 
<LI{e\u^-\s, y)-«("-l)(s, y)|2 


E|w('^)(s, z)-u^^-^\s, z)p 
r \ 1/2 


1/2 


< Ll\^Dn-iis,y)Dn-iis, z) 

< LlHr,_,{s). 
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Hence we have for 7 < a using Lemma 2.7 
Dn{t, x) 

<X^Llf Hn-i{s) f f Gt-six - y)Gt-s{x - z)f{y,z)dydzds 
Jo Jr‘‘ Jr<‘ 


< c.X^Ll 


/•* gn-l(s) 

/o 


ds. 


We therefore have 


Hn{t) < CiX^Ll [ 

Jo 


Hn-l{s) 

{t - 


ds. 


We now note that the integral appearing on the right hand side of the above 

display is finite when d < a/fJ. Hence, by Lemma 3.3 in Walsh [32], the series 
1 

J2'^=o converges uniformly on [0, T], Therefore, the sequence {u„} con¬ 

verges in and uniformly on [0, T] X and the limit satisfies (1.7). We can 
prove uniqueness in a similar way. 

We now turn to the proof of the exponential bound. Set 

A{t) := sup E|ut(a;)p. 

We claim that there exists constants C 4 , C 5 such that for all t > 0, we have 

A{t) < C 4 + C 5 {XL^f ds. 

The renewal inequality in Proposition 2.12 with p = {a — ^(3)1 a then proves 
the exponential upper bound. To prove this claim, we start with the mild 
formulation given by (1-7), then take the second moment to obtain the following 


E|ut(a;)p = \{gu)t{x)\^ 


fl 

Jo Jr" 


+ X 

h Jr'^xR'^ 
= h+h. 


Gt-six, y)Gt-six, z)fiy, z)E[aiusiy))a{usiz))]dydzds 


(4.1) 


Since uo is bounded, we have Ji < C 4 . Next we use the assumption on a 
together with Holder’s inequality to see that 


E[cr(u42/))cr(-«s(2;))] < LlE[usiy)usiz)] 

< [E\usiy)\Y^^[E\usiz)\^]^/'^ 

< sup E|Ms(a;)|^. 


(4.2) 


Therefore, using Lemma 2.7 the second term I 2 is thus bounded as follows. 


I 2 < 


coiXL^f f 

Jo 


{t - 


ds. 
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Combining the above estimates, we obtain the required result in the claim. 


□ 


4.2 Proof of lower bound in Theorem 1.10. 

The starting point of the proof of the lower bound hinges on the following 
recursive argument. 

= \{gu),{x)\^ + x^ f [ 

Jo J«.‘‘xR'^ 

G(t - Si, X, zi)G{t - Si, X, z[yE,[a{usi{zi))a{us^{z[)) f {zi, z[)]Azidz[Asi. 
We now use the assumption that (j{x) > la\x\ for all x to reduce the above to 
E|ut(a::)p 

>\iguUx)f+X^ll f [ 

Jo 

G{t - Si, X, Zi)G{t - Si, X, z[)E\usi{zi)usi{z[)\f{zi,z[)dzidz[dsi. 

We now replace the t above by t-\-t and use a substitution to reduce the above 
to 


>\igu\^i{x)f + xX f f 

Jo JR'^xR'^ 

G{t-Si, X, Zi)G{t-si, X, z[)E\uij^^^{zi)uij^^^{z[)\f{zi,z[)dzidz[dsi. 
We also have 


E\ui^^^{zi)ui^^^{z[)\ > \igu)i^^^{zi){gu)i^^^izi) \ + X'^ll [ [ 

Jo JR’^xR'^ 

G(si - S2, Zi, Z2 )G{si - S2, zi, Z2)E|'«t+s2 (z2)Wt+s2 (4) 1/(^2 , 4 )dz 2 dZ 2 dS 2 . 
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The above two inequalities thus give us 




Jo JR“ 


Jo JR'^xR'^ 

G{t-Si, X, zi)G{t-si, X, z[)E\u^_^_^^{zl)ui+s^{z[)\f{zl,z[)dzldz[dsl 

>\{gu)i^,{x)\^+XHl f [ 

Jo JR'^xR'^ 

G{t-Si, X, zi)G{t-si, X, z[)f{zl,z[){gu)^_^_^^{zl){gu)i_^_^^{z[)dzldz[dsl 

n pt + Si P 

G{t - si, X, zi)G{t - si, X, z[)f{zi,z[) / / 

;‘*xR'^ Jo Jr‘‘xR‘‘ 

G{si -S 2 , Zi, Z 2 )G{si - S 2 , Zi,' Z 2 )E\ui^^^{z 2 )ui^^^{z 2 )\f{z 2 ,Z 2 )dz 2 dz 2 ds 2 dzidz[dsi. 

(4.3) 

We set zq = z'q := X and sq ■= t and continue the recursion as above to 
obtain 

> \{gu)i^t{x)\^ 

+ f [ r [ ••• 

Jo Jr'^xr^Jo Jr^xR'^ Jo Jr^xr^ 

k 

]^G(sj_i - Si,Zi-i, Zi)G{si-i - Si,zl_^, zOf{zi,zl)dzidzldsi. 


2=1 

Therefore, 


(4.4) 


> \{gu)i+t{x)\‘^ 




k^l 


■ f [ \{gu)i+,jzk){gu)i+,^{z'i,)\ 

Jo JR’^xR'^ 


Jo Jr‘‘xR’^Jo Jr^xR"^ do 

]^G(si_i - Si,Zi-i, Zi)G{si-i - Si,z'_i, Zi)f{zi,zi)dzidZidsi. 


Proposition 4.1. There exists a tg > 0 such that for t > to, 

°° / i \ k{a-10)/a 

E|u(t + to: a;)|^ > ^ (A^l^ci) J whenever x € B{0,t^^°‘), 

where ci is a positive constant. 

Proof. We will look at the following term which comes from the recursive rela- 
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tion described above, 




k^l 


0 JO 


\{Qu)i+sJZk){Qu)lj^,^{z'k)\ 


■^xR'^ Jo JR'^xI 

]^G'(si_i — Si,2i_i, Zi)G{si-i — Si,zl_i, zl)fizi,zl)dzidz[dsi. 
i=l 

We can bound the above term by 




\ 2 k 


k 

Gi^Si—x Si^Zi—\^ Zi)G(^Si—-]_ Si^ Zj^_-^^ Zj^f (^Zi^ z^^dzidzj^dsi- 

i=l 

We now make a substitution and reduce the temporal region of integration to 
write 




\ 2 k 


ntik 


rt/k /. 

Jo v/ro 


ptjk 


^0 Jr^xr<^Jo ^r^^xr^^ Jo -'S(o,t^/“)xS(o,i^/‘^) 
k 

W_G{si,z^-l, Zi)G'(si,z-_i, z-)f{zi,z-)dzidzidsi. 

i=l 

We will further reduce the domain of integration so the function 

k 

Y[G{si,Zi-i, z,)G(si,z'_i, zl)f{zi,zl), 

i=l 

has the required lower bound. For z = 0, • • • , fc, we set 

Zi G B{x, sf^“/2) n B{zi-i, sf^“) 


and 


z'i G B{x, sf/“/2) n B{z'i_i, 


P/oi\ 


We therefore have fzi — z'j < \zi — Zi-\\ < and \z[ —z[_t\ < We 

use the lower bound on the heat kernel to find that 


„/5/a 


„/5/a 


Y[G{si,Zi-i, Zi)G{si,zl_^, zl)f{zi,zl) 


— k'j^jc 


^ k -€ 

c TT 1 


n 




2 ^dfo 


for some c > 0. We set A := B{x, s^^“/2)ni3(zi_i, sf^“) and^' := B{x, Si^°‘/2)n 
B{z[_i, We will further choose that sf^“ < and note that \Ai\ > 

and l^'l > We therefore have 
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OO ~t/fc . . 

/ / / 

JO JR'^xR'^JO JR 


f‘tlk 


R'ixR'^ Jo JB(0,t'3/“)xB(0.t'3/°) 
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We now take time large enough and use Lemma 2.11 to complete the proof 
of theorem. □ 


4.3 Proof of Theorem 1.12. 

The proof of this theorem is exactly as that of Theorem 1.8 and it is omitted. □ 


4.4 Proof of Theorem 1.13. 

Proof. We will make use of the Kolmogorov’s continuity theorem. Therefore we 
consider the increment E|itt+;i(a::) — Ut{x)\P for h G (0,1) and p>2. We have 


ut+h{x) - ut{x) = / [Gt+/i(a; - y) - Gt(a: - y)]uo(2/)dy 
Jr<^ 

a / [ [Gt+h-six - y) - Gt-six - y)]a{usiy))F{dsdy) 

Jo JR<i 


+ • 


+ A 




Gt+h-six - y)aiusiy))F{dsdy). 


(4.5) 


The first term, Gt(x — y)uo(y)dy is smooth for t > 0. This essentially 
follows from the fact that under the condition on the initial condition, we can 
interchange integral and derivatives. We will therefore look at higher moments 
of the remaining terms. Recall that we can use the similar ideas in the proof of 
Theorem 1.10 to show that sup 3 ,gRd E|ttt(x)|/’ is finite for all t > 0. We use the 
Burkholder’s inequality together with Proposition 2.6 to write 


28 
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(4.6) 


< cih 


p{l-0^/a) 


Similarly we have 


p 


E| / / Gt+h-Ax - y)(^iusiy))F{dsdy)\P 


t JR'’- 



<C2\ / / Gt+h-six-y)Gt+h-six-z)f{y,z)dsdydz\P/^ 


(4.7) 


p+h !< 1 

<C2\ / [Ep{-,z\^nt + h-s)^)f-^d^ds\P/^ 

Jt Jr'^ It I ' 



Combine the above estimates, we see that 


E\ut+h{x) - ut{x)\ 



Now an application of Kolmogorov’s continuity theorem as in [4] completes the 


proof. 


□ 
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